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We studied the non-equilibrium short-time dynamics of a spin- 1/2 chain with Dzyaloshinskii- 
Moriya interactions after a sudden quench by a transverse field. We found that inhomogeneous 
spin-spirals with opposite chiralities propagate from the edges towards the center of the chain. This 
propagation is accompanied by quantum spin oscillations which decay asymptotically with time. A 
theoretical description of this phenomenon is given to a good accuracy with the help of numerical 
calculations with the adaptive time-dependent density matrix renormalization group algorithm. 



I. INTRODUCTION 

In recent years, many attempts have been made to 
understand the non-equilibrium dynamics of closed in- 
teracting quantum systems^. These efforts are mainly 
motivated by the high degree of isolation and control of 
microscopic parameters which can be achieved experi- 
mentally in systems such as ultracold atoms trapped in 
optical lattices 2 -, which are currently being used to "en- 
gineer" quantum spin systems with desired properties^ — 
with the aim of simulating the behavior of emblematic 
systems in quantum statistical physics and discovering 
potential technological applications. Among the most 
studied problems which have been addressed, the dy- 
namics of low-dimensional systems following a quantum 
quench has received special attention. This is one of the 
simplest ways to drive a system far from equilibrium. It 
consists of preparing the system in the ground state of an 
initial Hamiltonian and then suddenly (or slowly) chang- 
ing the interactions at a given time, by varying an experi- 
mentally tunable parameter, with the result of letting the 
system evolve unitarily under a final Hamiltonian. As 
an example, one might mention the quenching dynam- 
ics by a transverse magnetic field (tunable parameter) 
in systems like the XY chain 6 9 or the quantum Ising 
chain-i^— , which can be simulated in ultracold atomic 
systems^. 

Quenching by a transverse field is a typical proto- 
col for studying non-equilibrium dynamics, since apart 
from being easy to implement in real situations, it al- 
lows the possibility of taking the system (in some cases) 
through a quantum phase transition, involving macro- 
scopic changes in the many-body state at the final point, 
as in the quantum Ising chain£. In more general quench- 
ing scenarios, the universal properties of the dynamics 
near quantum critical points as well as in generic gapped 
and gapless systems have been investigated thoroughly^; 
in some cases it is possible to make universal predictions 
of scaling laws for relevant physical quantities. Also, 
the question of whether the system relaxes to a quasi- 
stationary state in the long-time limit after the quench 
and, if it does, how to characterize its physical properties, 
has been addressed with a great interest. In particular, 
the concept of thermalization of a quantum many-body 



system following a quench has acquired importance; it 
has been conjectured that the asymptotic state can be 
described by a generalized Gibbs ensembl o 15 i 16 if the sys- 
tem is integrable, and by some effective thermal distribu- 
tion in the case of a generic system. Another main line 
of research, which we follow in this paper, is to study 
the characteristics of the time evolution just after the 
quench (short-time dynamics). This is of special inter- 
est, for example, in quantum computation scenarios^ in 
which gate operations require a precise knowledge of the 
immediate dynamical response of the system. Experi- 
mental examples following this line of research are the 
oscillations^ or the dephasingi^ of the matter wave field 
of a Bose-Einstein condensate, features which are also 
observed in the antiferromagnetic order parameter of a 
system of ultracold atoms in an optical lattice with su- 
perexchange interactions^. 

There are many spin systems which serve as proto- 
types for studying the short-time non-equilibrium dy- 
namics following a quantum quench. In this paper we 
study a spin system with exchange-relativistic interac- 
tions which, according to our knowledge, has not been 
investigated so far in this context. These kind of interac- 
tions arise when Anderson's superexchange mechanism 
is extended to include spin-orbit coupling. As is well 
known, the spin-orbit interaction (which couples the spin 
to the crystal lattice, giving rise to easy and hard mag- 
netization axes) is gaining considerable attention due to 
the possibility of manipulating spins in solid state sys- 
tems solely by electric fields 20 and/or perhaps by applied 
local strai n 21 ! 22 . In low-dimensional systems, the lack of 
structural inversion symmetry which is induced by sur- 
faces and interfaces, combined with the spin-orbit cou- 
pling leads to anisotropic antisymmetric superexchange 
or Dzyaloshinskii-Moriya interaction (DMI) 23 i 24 , which 
is responsible for the presence of weak ferromagnetism 
in a variety of antiferromagnetic compounds^. One of 
the most striking consequences of this spin interaction is 
the appearance of homochiral spin spiral magnetic struc- 
tures, as have been observed recently along fixed cryst al- 
lographs directions (easy-axis directions in some cases) 
in ferromagnetic and antiferromagnetic ultrathin films 
using spin-polarized scanning tunneling microscopy^ —. 
These one-dimensional spin structures have been postu- 
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lated to play a crucial role in spintronic devices, since 
spin-polarized currents flowing through them will exert 
a spin- tor que on the chiral magnetic structure, causing 
a variety of controllable excitations 2 ^. They are also the 
most likely candidates to host ferroelectricity, according 
to recent experimental results^ - — . In the context of 
ultracold atoms in optical lattices, it remains to be dis- 
covered whether the DMI can be artificially implemented 
and controlled, as can be done with the symmetric su- 
perexchange interaction^. Due to the great interest in 
the consequences of this antisymmetric interaction, this 
is expected to be achieved soon, since the spin-orbit cou- 
pling has already been demonstrated recently in these 
systems of atoms 31 . 

Motivated by the above-mentioned potential applica- 
tions as well as by the search for fundamental knowledge, 
in this paper we study the short-time non-equilibrium 
dynamics in the simplest system with DMI: a spin- 1/2 
chain, after a quantum quench is carried out by means 
of a transverse field. To be more specific, we consider 
a quantum quench which consists of preparing the ini- 
tial state |^o) =|— » • • • —>) by applying a strong mag- 
netic field along the x-direction, and then turning it off 
abruptly at t = 0. The subsequent unitary evolution, 
\ipt) — exp(— iHt)\ipo)i is governed by the final Hamilto- 
nian (we set h = 1 throughout this paper) 

H = ^[D{sjsy +1 -sysj +1 ) 

(i) 

+ J (SjSj+i + SjSj +1 + A SjSj +1 ) ], 

where the usual DMI in the z-direction is taken together 
with the Heisenberg-type symmetric exchange coupling 
(EC) with anisotropy A (also known as the XXZ ex- 
change interaction). As an example of a real quasi-one- 
dimensional compound described by this Hamiltonian we 
mention KCUF3, which is an antiferromagnetic spin- 1/2 
system with DMI and almost isotropic EC—. The crit- 
ical behavior of this system will be briefly discussed in 
section ITTT1 Meanwhile, we simply mention that the DMI 
makes the phase diagram rich, with two critical points 
(as in the XXZ model) dependent on D, and restores 
the otherwise spoiled ground-state entanglement due to 
anisotropy (which tends to align the spins), by means 
of enhanced quantum fluctuations 3 - 3 -. With respect to 
some dynamical and thermal features, this model has at- 
tracted attention in the quantum information commu- 
nity; for instance, the striking phenomenon of finite-time 
disentanglement— has been predicted for the two-qubit 
system prepared in the Werner stated, and teleportation 
of states has been studied using the model in thermal 
equilibrium as a quantum channel 36 . 

In most of the numerical calculations in this paper, we 
use the adaptive time-dependent density matrix renor- 
malization group algorithm (t-DMRG) 37 i 38 , which has 
proved to be a powerful method for studying the dy- 
namics of one-dimensional quantum many-body systems 
with short-ranged interactions. As a general result, we 



show that the DMI gives rise to inhomogeneous spin- 
spirals with opposite chiralities, which propagate from 
the edges towards the center of the chain until they meet 
each other, and then complex interference patterns arise. 
The propagation is accompanied by quantum spin oscil- 
lations that decay asymptotically with time. It is the aim 
of this work to theoretically describe these oscillations, as 
well as the nature of the propagation of the spin-spirals. 
The paper is organized as follows: in section [III we study 
in detail the role played only by the DMI in the for- 
mation of these spirals. For this, we use results from the 
two-spin quantum dynamics and many-body classical dy- 
namics to interpret the observations from the quantum 
model. In section IIII1 we add the exchange interactions 
and investigate how the results from section [II] are mod- 
ified. In section HVl we make some comments about the 
long-time dynamics in our system; and finally in section 
IVl we give the conclusions. 

II. EFFECTS OF THE DMI 

Although the DMI is usually expected to be negligible 
(due to its relativistic nature) compared to the EC, there 
are instance o 25 i 26 where it is strong enough (D / J ~ 0.25) 
to create remarkable spin-spiral structures. In this sec- 
tion, we take the limit J = in order to investigate the 
effects on the dynamics due only to the DMI. Then Eq. 
(PQ) is reduced to 

n ™ = D J2(S*S« +1 -S]S* +1 ). (2) 

3 

This model is exactly solvable, since it can be mapped 
onto the XX model: H xx = DJ2j (S?S? +1 + S] +1 ) 
(which is well known to be equivalent to a set of free 
spinless fermions) by rotating the basis at each site j by 
an angle j\ around the z-axis. Moreover, its eigenstates 
are non-equilibrium or spin-current-carrying stationary 
states of the XX model at zero temperature, since the 
macroscopic spin-current J M = Yj Jj 4 (quantized along 
the direction of D) defined through the continuity equa- 
tion 

S* = i [H™,Sf\=J? _jM^ (3) 

can be found 3 -^ to give J M = H DM and, given that this 
spin-current is a conserved quantity in the XX model, it 
follows that [H XX ,H BM ] =0. The time evolution result- 
ing from several initial states in the XX model such as 
the domain- wall state at the center of the chain22r— , the 
Neel state and the spin-density wave stated 2 - has already 
been solved analytically, yielding a variety of phenom- 
ena such as quantum fronts in which quantized steps of 
spin flips are transported and oscillatory dynamics of the 
order parameter. Our initial spin-polarized state in the 
x-direction |^ > oc (||)+ ||)) ® (|t>+ ID) • • • ® (|t>+ ID) 
contains eigenstates of the XX model such as the fully 
polarized states |tt * ' * t) and • • • |), and also many 
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FIG. 1. Time evolution of the components of magnetization 
for several values of the t-DMRG truncation parameter m in 
the arbitrary case with DMI and EC for which D/J = 0.48. 
(a) x-component. (b) ^-component, (c) ^/-component. 



excitations with a large number of domain walls (the 
Neel state being one of those); thus it is a highly non- 
equilibrium state at zero temperature. When this initial 
state is evolved under the XX Hamiltonian, there is no 
signature of dynamical formation of spin-spirals. Oth- 
erwise, when it is evolved under the DM Hamiltonian, 
the trajectory in the Hilbert space of the resulting state 
is different, and the formation of the spin-spirals is ob- 
served. Since there are no driving fields to excite this 
wave phenomenon, the open boundary conditions in our 
finite system play a crucial role in the establishment of 
these spin-currents. In particular, since (^o| J M \ipo) = 0, 
we expect to encounter a positive spin-current propagat- 
ing to the right from one boundary, and a negative spin- 
current propagating to the left from the other boundary. 
Although in principle the real-time evolution can be in- 
vestigated analytically in this case, it is a highly non- 
trivial problem; therefore we use the t-DMRG algorithm 
to inquire further into this. 



A. Many-body quantum dynamics 

We have used the t-DMRG algorithm for a chain with 
N = 200 sites, with a Trotter time step St = 0.05, and 
keeping m = 40 — 100 density matrix eigenstates to rep- 
resent the truncated sector of the evolved wave function. 
To test the convenience of this choice of simulation pa- 
rameters, we show in Fig. [I] the time evolution of the 
magnetization M a = (l/N) ^ZjiS^) (a = x, y, z) for sev- 
eral values of m, in an arbitrary case with J ^ in which 
strong long-ranged correlations are expected to lower the 
accuracy of the algorithm compared to the case J = 0. 
As can be seen, for values m > 60, the m-dependent error 
in the chosen time window is of order 10 ~ 7 in the worst 
case, which is a satisfactory error for DMRG users. In 
the following calculations, we used m = 60 (unless oth- 



FIG. 2. Snapshot of the local magnetization and nearest- 
neighbor spin-spin correlations in the x-direction at the arbi- 
trary time Dt — 80. (a) As can be seen, inhomogeneous spin- 
spirals with opposite chiralities propagate, in the xz-plane, 
from the edges towards the center of the chain, accompanied 
by spin oscillations. A classical picture of the left-rotating 
spiral arising from the left end is drawn when Dt = 2 (only 
the directions of spins are taken into account), (b) The spirals 
propagate with the same velocity with which nearest-neighbor 
correlations move. 



erwise stated) with the aim of saving computer resources 
and speeding them up. Note that when the value of J is 
not shown in the figures, it is assumed to be zero. 

Returning to the case of relevance for this section 
(J = 0), we show in Fig. [2] a "snapshot" of the local mag- 
netization along the chain as well as the nearest-neighbor 
spin-spin correlations CJj + i = (SjSj +1 ) — (5J)(5J +1 ) in 
the direction of the polarization of the initial state, at 
the arbitrary time Dt = 80. As can be seen, inhomoge- 
neous spin-spirals with opposite chiralities propagate (in 
the xz-plane) from the edges towards the center of the 
chain, as expected from the spin-current considerations 
given above. This propagation is accompanied by spin os- 
cillations around the z-axis. A sketch of the left-rotating 
spiral propagating from the left end of the chain is drawn 
when Dt = 2. We verified that the chirality (sense of ro- 
tation) of the spin-spirals reverses when the sign of D 
is reversed, which is an inherent property of the DMI2 7 -. 
We also followed the wave fronts, and found a constant 
velocity of propagation v = 1 (in units of aD, a being 
the lattice constant). As can be seen in Fig. [2fb), this 
is the same velocity with which nearest-neighbor corre- 
lations propagate from the edges, in accordance with the 
light-cone propagation which is expected for spin models 
with sufficiently local interaction s 43 ! 44 . 

It is also interesting to investigate the time evolution 
of long-ranged correlations, which we calculate as is usu- 
ally done in the DMRG framework 4 ^ C x r = (SfS?) - 
(Sf >(5?>, where i = iV/2 + l- r and j = N/2 + r (when r 
is increased, i and j move outwards symmetrically from 
the center towards the edges of the chain). We found 



that long-ranged correlations are quite small in this case: 
they are always below O(10 -3 ) during the spiral propa- 
gation, compared to the nearest-neighbor ones which, as 
may be seen in Fig. Efb), are O(10 _1 ); this is why they 
are not shown. We shall see later that when J ^ 0, strong 
long-ranged correlations develop throughout the system, 
making the dynamics more intricate. Meanwhile, let us 
first try to understand the nature of the results obtained 
so far. 



B. Two-spin quantum dynamics 

We begin by considering the quantum dynamics of a 
two-spin system, which should contain relevant informa- 
tion about the behavior of a pair of nearest-neighbor sites 
in the many-body system. In order to simplify the de- 
scription of states, we rotate the basis temporarily by 
changing x o z in Eq. (pQ), which leads to a DMI in 
the x-direction and a fully polarized initial state in the 
z-direction. Note that this operation does not change ex- 
pected values of observables. It is also convenient to in- 
troduce the Bell states in this new representation, defined 
as \t/s) = (|U)± |lt))/V2 and \t/s) = (|ft)± \U))/V2. 
Then we show in Fig. [3fa) the energy levels and the sta- 
tionary states of the two-spin model, with ujr = D/2. 
Energetically speaking, this is equivalent to a two-level 
magnetic system with a spinless intermediate state. In 
the present case (J = 0), when the initial state |^o) =|tt) 
is prepared, the evolved two-spin state reads 
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\t)+cos{uj R t)\s) -sm(u R t)\s)}, (4) 



which shows Rabi oscillations between the states |tt) an d 
i.e., the ^-components of local magnetization for 
the spin pair oscillate in phase (5f ) = = \ cos(u R t) 
with the Rabi frequency ujr. Remarkably, the most prob- 
able state to find after an energy measurement of the 
spin pair (during a cycle) is the spinless intermediate Bell 
state \i), so, at any time, spin measurements in the di- 
rection perpendicular to D will disperse around zero, as 
indeed is seen in Fig. EJa). 

Returning to the original representation, we have a 
two-spin system with DMI in the z-direction and initially 
prepared in the state |^o) =|— The time evolution 
shows Rabi oscillations between the states !—>•—>•) and 
with (Sf) = (S%) = \ cos(ojRt). A natural ques- 
tion then arises: how can we generalize these results to 
many sites? It is reasonable to assume that, at least for 
times Dt < 1, the evolved state is close to 
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FIG. 3. Energy levels and stationary states of the two-spin 
model with DMI. (a) In the absence of EC, the first excited 
state is two- fold degenerate, (b) In the presence of EC, the 
first excited state remains degenerate, but the energy gaps 
are modified. 



end of the chain (negative spin-current). This assump- 
tion is motivated by the fact that for Dt < 1, spin-spin 
correlations have not crossed the boundaries of the two- 
site blocks at the ends of the chain, so the state can be 
thought of as factorized into "non-interacting dimers": 
| • • | | • • | • • • | • • | | • • | . i.e., a product of two- 
spin solutions. The problem with this assumption is that 
it does not take into account the spiral propagation; i.e., 
since the perturbations originated at the chain bound- 
aries propagate with a finite velocity, it is clear that the 
spins within these dimers must oscillate with different 
phases, so the state in Eq. ([5]) should be modified to 
accomplish this. 

A convenient way to account for the spiral formation 
is to decouple the two-spin dynamics by introducing the 
reduced density matrix pf = Tr2|'0^ = )2s('0t = | for the spins 
that have been reached by the perturbations, and po = \—> 
)(— H for those that have not. Then, an easy calculation 
leads to pf = \X + cos(ujRt)S x ± | sin(Dt)S z and po = 
\l> + S x , where X is the 2x2 identity matrix. Note that 
this decoupling of the two-spin dynamics preserves the 
relation (S x ) = Tr(pfS x ) = ±cos(u> R t). With this, we 
expect that the evolved state will be close to 
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3 v > vjy j 

j=N/2+l 



(6) 



where A^fi) = [1 — H(rj)]po + H(rj)p^. selects the ap- 
proximate state of the site j at time t, with H being the 
>) \ipt~)2s, (5) Heaviside step function and Tj = t—(j — l)/D (tj = the 



where \^f)2s is the state in Eq. (j4]) written in the original 
representation and \^)2s is the corresponding state ob- 
tained when changing D <B> —D, in order to account for 
the opposite chirality of the spiral arising from the right 



moment at which the perturbation hits the site j, mea- 
sured from the chain boundaries). We tested the validity 
of the ansatz in Eq. ([6]) by calculating the exact time evo- 
lution of local magnetization for a system with 14 sites, 
as shown in Fig. 2J As can be seen, the ansatz works 
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FIG. 4. Snapshots of the local magnetization for a system 
with 14 sites, calculated exactly: (a), (c); and with the ansatz 
in Eq. JB}: (b), (d), for two different times. The black- filled 
circles are values of (Sj), the gray- filled ones are those of (Sj), 
and the white- filled squares are values of (Sj), as in Fig. [2ja). 
Solid lines are guides for the eyes. 




Dt 



FIG. 5. Time evolution of the x-component of local mag- 
netization for the first spin in a system with 14 sites, (a) 
Comparison between exact results with those obtained from 
the ansatz in Eq. ©; solid and dashed lines are guides for 
the eyes, (b) Finer curve obtained with t-DMRG for times 
Dt<2. 



very well for times Dt < 1 as expected from the argu- 
ments motivating Eq. (|5]) . A closer look at the single-site 
dynamics is shown in Fig. [5j up to the time just before 
the propagating spirals interfere with each other at the 
center. Fig. [5jb) shows the comparison between the 
x-component of the local magnetization calculated with 
the ansatz and with t-DMRG using a Trotter time step 
St = 0.01 and m = 100 (it would take too much time to 
sweep the exact curve). As is evident, the ansatz gives 
satisfactory results even for times Dt < 2. For larger 
times, the assumption leading to Eq. (j5j) is no longer 
valid, since the perturbations have reached the third sites 
(measured from the chain boundaries), and three-body 
effects then arise. In general, as nearest-neighbor spin- 
spin correlations propagate towards the center, many- 
body effects become important and the single-site dy- 
namics deviates remarkably from what is expected from 
Eq. (|6j). As we shall discuss next in more detail, the most 
remarkable many-body effect is the onset of an "emer- 
gent" Rabi frequency of oscillations uj^ in the single-site 
dynamics; these oscillations decay with time as can be 
anticipated from Fig. [5ja). 



C. Many-body classical dynamics 

The simplest method to investigate many-body effects 
in our system is by appealing to classical considera- 
tions. At the limit of commuting dynamical variables, 
the Heisenberg equations of motion for the spin compo- 
nents turn into the Landau-Lifshitz equations (without 
damping) for the dynamics of classical spin vectors Sj 



precessing around effective fields Bj \ 



dSi 

dt 



-Sj x Bj, with Bj 



on 



DM 



dSi 



(7) 



In polar coordinates Sj = \ cos 9j , where 9j is the prop- 
agating spiral profile, and if we assume 5j = 0, as is 
evident from Figs. [2ja) and |4j we obtain the system of 
differential equations 



0j = (jjf>(cos9j+i — cos0j-i) 



(8) 



Let us consider the classical two-spin model first. Just 
after t = 0, the first two spins are almost collinear and 
thus 0i « #2- With this approximation, we have 0\ = 
cjrcos#i and 62 = —ur cos 62, with the solution 



25f it) = 2Sf (t) = sech(c^). 



(9) 



This is plotted in Fig. [6fa) with a dash-dotted line. It 
suggests that the classical spin pair (decoupled from the 
many-body system) moves in phase, as is the case in the 
quantum two-spin solution. However, as may be seen 
from the figure, the classical dynamics of the spin pair do 
not display any kind of oscillation. As we shall see next, 
this is not the case in the many-body classical system. 

The first step for solving the many-body classical dy- 
namics is to numerically solve the system of differential 
equations in Eq. (j5J). From this system we observe that, 

since 0\ = uor cos 62 and On = — cjr cos0/v-i> ^ ne or i" 
gin of spin-spirals with opposite chiralities at the chain 
boundaries is a straightforward consequence of the clas- 
sical open boundary conditions (see the change of signs). 
The solution of these equations for the first two spins 
along the chain is shown in Fig. [6] with thick solid lines. 
As can be seen, the classical solution is able to predict 
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FIG. 6. Time evolution of the x-component of local magneti- 
zation for (a) the first spin and (b) the second spin along the 
chain. The black- filled circles are the results from t-DMRG. 
The dash-dotted line in (a) is the result of the classical two- 
spin model described by Eq. (|9]). The solid line is the numeri- 
cal solution from the many-body classical dynamics governed 
by Eq. ©, and the dashed line is the prediction from the 
ansatz in Eq. ©. 



the decaying oscillations over time, although the decay 
is far slower than the data from the quantum model (re- 
sults from t-DMRG). Also, the observed frequency of the 
oscillations uj^ is very well reproduced by the classical 
solution. Comparing that with uor obtained from the 
ansatz in Eq. ((6]) (dashed lines), which neglects many- 
body correlations, we have uo^ = 2ur. Thus we conclude 
that the Rabi oscillations in the quantum two-spin sys- 
tem persist in the many-body system, but with a differ- 
ent "emergent" frequency. This kind of behavior has ap- 
peared recently^ 2 - in the context of quantum quenches in 
the anisotropic spin-1/2 Heisenberg chain. There, when 
the initial Neel state is prepared and subsequently 
evolved under the two-spin XX Hamiltonian, Rabi oscil- 
lations (between the states |ti) and ||t)) are observed 
in the staggered magnetization (antiferromagnetic order 
parameter) with a frequency ujr = J. However, when 
a sudden quench is performed from the many-body Neel 
state to the critical phase in the XX model, the order 
parameter still displays Rabi oscillations, but with the 
emergent frequency = 2cjr, as in the present case. 

The classical single-site spin oscillations coming from 
the many-body numerical solution take place with the 
same frequency as in the quantum model, although a 
slight dephasing between the two results can be observed 
in Fig. [6j This is because the classical solution is not 
able to reproduce the velocity v = aD with which the 
spin-spirals propagate (see e.g. the beginning of the cor- 
responding curves in Fig. EJ^b)). We said before that the 
present quantum model is equivalent to a set of free spin- 
less fermions. In fact, the resulting dispersion relations 
are the same as those of the XX model, but shifted by 



the wave number k a = x , since 

h bm = sin (fc a ) c t Cjb _ ^^ e - ia(fe -H (10) 

k k 

where c\ (c&) are creation (anihilation) operators of free 
spinless fermions with momentum k. Within the ther- 
modynamic limit, the second term of this expression van- 
ishes, and the group velocity of propagating excitations is 
given by Vk = aDcos(ak). Then, the propagating spin- 
spirals can be associated with free-fermion excitations 
with zero total momentum (zero energy modes) giving 
vo = v = aD. Hence, since this velocity comes from a 
purely quantum phenomenon, it is reasonable that clas- 
sical considerations do not reproduce it. 

In summary, we have learned from this section that 
the DMI term only, being exactly spin-current (quan- 
tized along the direction of D), drives the formation of 
spin-spirals which can be described qualitatively in classi- 
cal terms, but needs input from quantum theory to fully 
account for the observed behavior. In particular, the 
two-body quantum dynamics are helpful for interpreting 
the more complex quantum many-body results. In the 
following section, we investigate how these results are 
modified when isotropic EC is added to the system. 

III. ADDING EXCHANGE INTERACTIONS 

In order to investigate the effects of EC on the non- 
equilibrium dynamics studied above, we consider some 
limits of the anisotropic Heisenberg exchange interaction 
in Eq. (pQ). First note that, by rotating the basis on odd 
(or even) sites by an angle tt around the z-axis, the an- 
tiferromagnetic case of the Hamiltonian (J > 0) can be 
mapped onto the ferromagnetic case (J < 0) with op- 
posite sign of the anisotropy A. In this paper we take 
J > as our overall energy scale (measuring time in 
units of J -1 ), and give a brief overview of the quantum 
phase diagram as the non-dimensional parameters A and 
D j J are varied. For numerical calculations we set J = 1. 
The critical behavior of the (equilibrium) ground state 
of the model in Eq. (pQ) is well known^Zr— . For A < 
— \/l + {D / J) 2 , the ground state of the system is ferro- 
magnetic. For -y/l + (D/J) 2 < A < y/l + (D/J) 2 , the 
energy spectrum is gapless and the ground state does 
not support any kind of long-ra nged order and is thus 
known as a spin fluid. For A > yl + (D/J) 2 , the ground 
state is antiferromagnetic (Neel ordered). In both the 
spin fluid and Neel phases, the helical spin structure in 
the £2/-plane is obtained, which is one of the character- 
istic properties of the DMI. At the limit A — » oo, the 
Ising model with DMI is realized. The critical proper- 
ties and entanglement of this model have been studied 
recent l y 5Q i 51 , revealing two quantum phases: the antifer- 
romagnetic one and the saturated chiral one, separated 
by the critical point D/J=l. At the limit A = 0, the 
XX model with DMI is obtained. Note that by rotating 
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FIG. 7. (a) Snapshot of the local magnetization at the ar- 
bitrary time Jt — 70, in the model of Eq. (pQ) with A = 1 
and D/J = 0.8. A classical picture of the left-rotating spiral 
arising from the left end of the chain is drawn when Jt = 2, 
projected on the xz-plane (only the directions of spins are 
taken into account), (b) Nearest-neighbor spin-spin correla- 
tions, (c) Spin-spin correlations measured symmetrically from 
the center of the chain. 



FIG. 8. Time evolution of local magnetization for the first 
spin along the chain for the arbitrary value D/J — 0.8. In- 
set: Period of oscillations of the x-component of local mag- 
netization as a function of D/J. The solid line is the curve 
Tr = 2tt/ujr, the dashed one is the curve T u — 2k/uj Ui the 
dash-dotted line is the curve Tu = 2ty/(oj u — ^d), and the 
symbols are the results from t-DMRG. 



the basis at each site j by an angle (j — 1) tan -1 (D/J) 
around the z-axis, it can be shown^ that the DMI term 
of the Hamiltonian can be removed, leading to a XX 
model with a renormalized exchange interaction energy 
\/J 2 + D 2 . In the isotropic case (A = 1), the DMI term 
can also be removed from the Hamiltonian^ (in fact, this 
can be done for any value of A; see e.g. [47]) after a simi- 
lar canonical transformation. The resulting Hamiltonian 
is that of the XXZ Heisenberg model 

h xxz = yj[ j**(sjsj +1 + §?s? +1 ) + rsfs z j+1 }, (ii) 

3 

with planar exchange energy 

jxx = v /j2 +jD 2 ? and 

anisotropy energy J z = J. Since J z / J xx < 1 for any 
non-zero value of D, the isotropic system with DMI will 
be in a spin-fluid phased, as can be inferred from the 
discussion above. 

The relaxation dynamics of the antiferromagnetic or- 
der parameter after several quantum quenches in the 
XXZ model has been investigated in detail recently^ 2 -, re- 
vealing a new mode of many-body dynamics in the effec- 
tively interacting case of finite anisotropy (J z / ' J xx in Eq. 
([TT]) ): a exponential decay of the staggered magnetization 
with (or without) oscillatory behavior. In a mean-field 
approximation, a rough picture of the dynamics can be 
obtained by introducing pseudo-spins in /c-space. Start- 
ing from the Neel state (equivalent to all pseudo-spins 
pointing in the x-direction), the pseudo-spins begin to 
precess around a (fc, t)-dependent Zeeman field (static, 
and in the z-direction in the effectively non-interacting 
case of the XX model). Since the field (and hence the 
precession frequencies) varies continuously from pseudo- 
spin to pseudo-spin over a bandwidth, they gradually 



dephase, causing the oscillations to decay with time. A 
similar quasi-classical picture (in real space) has been de- 
scribed recently in the study of coherent transport of spin 
currents in the spin- 1/2 Heisenberg chain in transversed 
magnetic fields^. Due to the connection with a real ma- 
terial like KCUF3, in the following discussion, we restrict 
ourselves to the isotropic case (A = 1) of the Hamilto- 
nian in Eq. (pQ) (equivalent to Eq. ([TT]) ). focusing on the 
description of the frequency of single-site spin oscillations 
by means of the two-spin quantum model. We verified 
that, at the Ising limit (A — >• oo), a similar method can 
also approximately describe the observed behavior . 

In Fig. we show a snapshot of the local magneti- 
zation and spin-spin correlations in the ^-direction at 
the arbitrary time Jt = 70 and D/J = 0.8. As can 
be seen, inhomogeneous spin-spirals with opposite chi- 
ralities propagate again from the boundaries towards the 
center of the chain. As shown in the classical picture 
projected on the xz-plane when Jt = 2, the velocity of 
propagation is, in general, greater than in the case with 
J = 0. As expected, the presence of EC promotes the on- 
set of strong long-ranged spin-spin correlations (see Fig. 
[7(c)). However, these correlations seems to be large for 
spin pairs (with respect to the center of the chain) that 
have not been touched by the spiral perturbations and 
diminish strikingly to zero for those that have undergone 
single-site oscillatory dynamics for a considerable time; 
that is, the propagating spin-spirals with opposite chi- 
ralities are almost uncorrelated. We also observe from 
Fig. [3(b) that the nearest-neighbor spin-spin correla- 
tions are strongly weakened compared to the case with 
J = 0, going to negligible values at (and near) the chain 
boundaries. A remarkable feature from Fig. [7(a) is that 
in the presence of EC, the ^/-component of local magne- 
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tization is no longer "quenched", as it was in the case 
with J = 0. A closer look at the single-site quantum 
dynamics is shown in Fig. [8j where the time evolution 
of local magnetization for the first site along the chain 
is calculated. We immediately see that beginning in the 
state |— >•), this spin precesses (or "spirals") clockwise to- 
wards the asymptotic state \\) with a frequency which 
we shall investigate next. We verified that the decay of 
oscillations is exponential. At the other boundary, the 
spin precesses counterclockwise towards the asymptotic 
state |^) with the same frequency. These perturbations 
are transmitted to the spins in the bulk apparently in 
a diffusive way, in contrast to the ballistic propagation 
which is observed in the case with J = 0. 

Based on the experience of the former section, we may 
want to consider the two-spin quantum dynamics of the 
model with DMI and isotropic EC. In Fig. [3fb) we show 
the energy diagram with ujd = J/2 and uj u = \/w 2 R + uj 2 d . 
In this case, the energy gap between the ground and first 
excitated state (which remains two-fold degenerate) is 
larger than in the case with J = 0. The x-component of 
local magnetization of the spin par oscillate in phase 

2(Sf) = cos(uJdt) cos(uj u t) + — sm(ujdt) sm(uj u t). (12) 

The nature of these oscillations depends strongly on the 

ratio uJd/^u = yl + (D / 'J) 2 which coincides with the 
anisotropy J z / J xx in Eq. (fTTj) . Then, when D/J ^> 1 
the leading behavior is governed by the first term of the 
sum in Eq. ([T2]) . This term describes fast oscillations 
with period T u = 2tt/uj u modulated by the cosine curve 
cos(cj^t), i.e, a beat-like pattern. According to Eq. ([IT]) , 
in this regime, the many-body system tends to be de- 
scribed by the XX model and we expect that these fast 
oscillations will resemble the Rabi oscillations found in 
the former section. To test these predictions, we cal- 
culate the period of oscillations (by averaging the time 
between sequential troughs) for the first spin along the 
chain as a function of D / J, using t-DMRG. This is shown 
in the inset of Fig. [51 There, we see how the dashed 
line (representing T u ) tends to coincide with the solid 
one (representing Tr = 2k/uur) in the regime D/J ^> 1, 
and both curves well describe the actual period of the 
single-site oscillations in the many-body system. On 
the other hand, when < D/J <C 1, the behavior in 
Eq. (fT2|) is approximately well described by the curve 
cos[(co> n — LJd)t] which displays oscillations with the long 
period = 2tt/(uj u — ujd)- This is represented in the 
inset of Fig. [8] by the dash-dotted line. At the limit 
D/J (isotropic Heisenberg model), the spin oscil- 
lations are not expected to be observed, as inferred by 
the divergence of the period of oscillations. In fact, our 
initial spin-polarized state in the x-direction is a station- 
ary state of this model and hence no time evolution is 
detected in local observables, as we have verified numeri- 
cally. It is curious to note that the period of spin oscilla- 
tions in the two regimes can be remembered easily if we 
make an analogy with a radiative electronic transition. It 



(a) 

(b) tttttttt„^ wul 

FIG. 9. (a) Initial spin-polarized state in the x-direction. (b) 
Approximate classical picture of the appearance of the state 
after the spin-spirals with opposite chiralities interfere at the 
center of the chain (see the text for deviations). 



is well known from elementary quantum mechanics that 
if an electron is undergoing a transition between station- 
ary states n and I, its position oscillates sinusoidally with 
an angular frequency (E n — Ei)/%. Thus, in the regime 
D/J ^> 1, we might think of the spins as "undergoing" 
radiative transitions between the ground and first excited 
state (see Fig. [3jb)) and hence oscillating with the pe- 
riod T< « 2ir/((jj u +ijJd), and in the regime < D/J <C 1, 
they might be thought of as undergoing radiative tran- 
sitions between the first and highest excited state and 
hence oscillating with the period T> w 2tt/(uj u — Ud). 

One may ask about the validity of using the quan- 
tum two-spin model to interpret results from the many- 
body system in the present case. In the former section 
we argued that in the absence of EC this can be done 
for a certain time interval, since long-ranged spin-spin 
correlations were almost "turned off" during the spiral 
propagation and only the transport of nearest-neighbor 
spin-spin correlations were responsible for the onset of 
quantum many-body effects, which make the predictions 
from the quantum two-spin model fail. In the presence 
of EC, we see that strong long-ranged spin-spin correla- 
tions develop in the system, and yet the single-site spin 
oscillations can be described approximately well by the 
quantum two-spin model. We believe that this happens 
due to the fact that only nearest-neighbor spin-spin cor- 
relations carry the information of relevant many-body 
effects and these correlations are very small in the pres- 
ence of EC, especially at the boundaries (see Fig. [3(b)). 
As noted from the inset of Fig. [5J this picture is very 
suggestive, since there are no "emergent" Rabi frequen- 
cies in the presence of EC; i.e., the frequencies predicted 
by the quantum two-spin model are those approximately 
observed in the single-site dynamics of the many-body 
system. This was not the case in the absence of EC, since 
nearest-neighbor spin-spin correlations were remarkably 
larger during the spiral propagation. 



IV. COMMENTS ON THE LONG-TIME 
DYNAMICS 

In the previous sections, we studied the short-time dy- 
namics of the spin-spiral propagation and oscillations in 
a spin chain with DMI and with(out) EC. By this, we al- 
ways meant the time evolution just after the quench and 
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before the two perturbations coming from the boundaries 
interfere at the center of the chain. Even though we are 
not going to go deeper into what happens after this, it 
is worthwhile saying something for the sake of complete- 
ness. As might be expected intuitively, when the spin- 
spirals with opposite chiralities encounter one another, 
complex interference patterns arise. We observed in all 
cases that the components of local magnetization trans- 
verse to the direction of the DMI take on very small val- 
ues as each propagating perturbation goes across the op- 
posite half of the chain from where it was created, which 
is reasonable, due to the "destructive" character of in- 
terference in the x?/-plane: spins in different halves of 
the chain have opposite direction of rotation (when they 
precess) in the x?/-plane before the propagation, and al- 
most uncorrelated spin-spirals reach the center. Then, 
beginning with our initial spin-polarized state in the re- 
direction as depicted in Fig. [9fa), we arrive after a long 
time at a state with an appearance approximately similar 
to that shown classically in Fig. [9fb), but with small fluc- 
tuations in the x?/-plane which support the propagation 
of residual spin waves. Since we could not access very 
long times with our numerical calculations, it was not 
possible to study properties in the very long-time limit. 

V. CONCLUSIONS 

We have investigated the non-equilibrium dynamics of 
a spin- 1/2 chain with Dzyaloshinskii-Moriya interactions 



in the presence (and absence) of Heisenberg exchange in- 
teractions after a quench by a transverse field. We found 
in all cases that inhomogeneous spin-spirals with oppo- 
site chiralities propagate from the boundaries towards 
the center of the chain until they meet at the center, 
producing complex interference patterns. This propaga- 
tion is accompanied by quantum spin oscillations that 
decay asymptotically over time. We showed that impor- 
tant information for the interpretation of the observed 
quantum behavior can be extracted from classical consid- 
erations of the spin system, as well as from the quantum 
two-spin model of interactions. Our results contribute 
to recent efforts to understand the characteristics of the 
non-equilibrium dynamics following a quantum quench in 
closed one-dimensional spin systems, and are expected to 
be helpful for possible potential applications where the 
Dzyaloshinskii-Moriya interaction plays a relevant role in 
low-dimensional systems. 
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